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Abstract 

In this paper, we explore the two-dimensional behavior of cellular automata with shuffle updates. 
As a test case, we consider the evacuation of a square room by pedestrians modeled by a cellular 
automaton model with a static floor field. Shuffle updates are characterized by a variable associated 
to each particle and called phase, that can be interpreted as the phase in the step cycle in the frame 
of pedestrian flows. Here we also introduce a dynamics for these phases, in order to modify the 
properties of the model. We investigate in particular the crossover between low- and high-density 
regimes that occurs when the density of pedestrians increases, the dependency of the outflow in the 
strength of the floor field, and the shape of the queue in front of the exit. Eventually we discuss the 
relevance of these results for pedestrians. 


I Introduction 

Modeling traffic systems has been a problem of growing interest in statistical physics during the last decade [U 

II [31 m. Indeed, these problems usually involve a large number of interacting driven agents that exhibit 
collective effects, and therefore belong to the realm of nonequilibrium statistical mechanics. Some examples of 
these collective effects are given by the spontaneous congestion of a highway [Slz] or the pattern formation 
observed in counter-propagating lanes of pedestrians [a El ED] or at intersections [81E1EI1E2I. Such effects can 
often be reproduced in minimal models inspired by statistical physics ESElEilESlEZlEHlEllin]- Many of 
these models belong to the class of cellular automata. In these models, the moving agents (cars, pedestrians, 
molecular motors, etc) are usually modeled as particles hopping on a lattice according to predefined rules, and 
interacting in particular through exclusion rules. 

Cellular automata are defined not only by the specification of the underlying lattice, and of the hopping rules, 
but also by the update order, i.e. the order in which the rules will be applied to the set of particles. It is known 
that update schemes can have a strong influence on the global dynamics of the system m- While random 
sequential update 0 is preferentially used in fundamental studies on out-of-equilibrium systems because of its 
closeness to continuous time dynamics, other updates with lower fluctuations are used in traffic applications. 
The most widely used is the parallel update [3], for which all particles are updated at the same time. In contrast 
with the random sequential update, the time step of parallel update can be given a physical meaning and be 
interpreted as a reaction time. 

Some other update schemes have been proposed, such as the random shuffle update [D2| and the frozen shuffle 
update [D3], both inspired by pedestrian applications. We shall emphasize that all these shuffle updates can 
be defined in terms of variables associated to the particles, and called phases. The phases determine the order 
in which particles are updated. Both the random and frozen shuffle updates were investigated in detail in one 
dimension [Ml US US nil 123 UHl US] • 

In this paper, we want to explore the properties of these updates when applied to a two-dimensional model. 
As a test case, we chose to consider the evacuation of a square room. Pedestrians are modeled by a simple 
so-called floor field model [m that we will describe more precisely in section 

*Now in Physics of Complex Systems, Weizmann Institute of Science, Rehovot 76100, Israel 
^In the random sequential update, one particle at a time step is chosen randomly and updated. 
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Figure 1: (a) Geometry of the room for L = 11, / = 5 and (b) possible hopping directions in two particular conhgurations. 
The dots represent particles and the grey cells are walls. The cells of the room have coordinates {x,y) with |x| < 1 and 
0 < 3 / < 21 + 1, and (0, 0) is the exit. The arrows represent the allowed hops for the particles in the centers. In the case 
of hybrid shuffle update which is defined in subsection |2.4[ the phase r stays the same for the hops with blue arrows and 
is redrawn in the case of red arrows. 


When the density increases, a crossover occurs between low- and high-density regimes. We will characterize 
this crossover, and the outflow in both regimes, including its dependency on the strength of the floor field. We 
will also discuss briefly the shape of the queue in front of the exit. 

Another aim of this paper is to illustrate by an example how it is possible to introduce a dynamics for the 
phases defining the update, and how this can modify the properties of the model. Eventually we shall discuss 
these results in view of pedestrian applications. 

This article is organized as follows. In section]^ we define the model in detail, and present the various shuffle 
update schemes used in the paper. In particular, we introduce some phase dynamics to define a hybrid shujfle 
update. Some theoretical and numerical results are presented in sections and where we study two different, 
low- and high-density regimes. In sectionj^ we turn to pedestrian applications and summarize some experimental 
results as well as modeling approaches from the literature. We discuss the most prominent features of evacuation 
flows that one would like to reproduce in models. In section we give the conclusion of this article. 

2 Model 

In this article we consider one of the simplest situations of pedestrians evacuating a room in the tradition of 
cellular automata. The room is divided into cells. Pedestrians are simply represented by ‘particles’ hopping 
stochastically from one cell to one of the neighbouring cells. Each cell can be occupied by at most one particle 
(‘simple exclusion’). Several evacuation models of this kind have been proposed [161 [TTl [30l [STJ [18] , in which the 
particles are usually updated in parallel, an update also traditionally used in car traffic problems. 

2.1 Lattice geometry 

We consider a square room oi L x L cells, where L = 2/ -|- 1 is an odd number, see figure for L = 11. The 
room has a unique exit located in the middle of a side, whose width is just the size of one cell. The cells of 
the room will be designated by their coordinates {x, y) with x = —I, + and y = 1,..., 2Z -|- 1, with a 

supplementary cell (0,0) representing the exit. 

In the following, particles will be allowed to hop from one cell to its empty neighbours, and a hop from (x, y) 
to (x', y') will sometimes be denoted as (x, y) —>■ (x', y'). This latter notation implicitly assumes that x, y, x' and 
y' are such that the hop is indeed allowed, i.e. {x,y) and {x',y') are neighbours, {x,y) is occupied and (x',y') is 
empty. 

2.2 Hopping probabilities 

At initial time and during the whole time evolution the simple exclusion constraint is verified, i.e. there can be 
at most one particle on each cell. 

We direct the particles towards the exit by using a static floor field [1^: The probability that each particle 
hops to a certain target cell is determined by the distance between the target cell and the exit of the room. 
We define a parameter k > 0 that quantifies the strength of the attraction towards the exit. We denote the 
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Euclideaijj norm by |r| = +y‘^ for r = {x,y). To each cell of the room we associate a weight 

u;(r) = e-'=l''l. (1) 

A particle on cell {x,y) chooses its target cell among all of its von Neumann neighbourhood (cells {x ± l,y) 
and [x, ?/ ± 1)) that are empty and the cell it already occupies. Among the possible targets, site r' is chosen with 
probability 

p(r') = Z~^w[r') (2) 

with the normalization Z = ?c(r) + where the summation runs over all the empty von Neumann 

neighbours of r. Figure (b) summarizes the possible hops of a particle. 

Finally, cell (0,0) is an exception to the rule (H): a particle standing on (0, 0) exits the system with probability 
1. In a larger simulation, its hopping probability would be determined by the downstream configuration. 


2.3 Shuffle update schemes 

When we use cellular automata, we need to specify an update scheme to fully define the model. In this paper, 
we shall exclusively consider shuffle updates. 

Shuffle updates are sequential updates, i.e. particles are updated one after the other. Here we chose to 
formulate these updates by associating to each particle i a phase Ti (0 < < 1). Then the initialization 

procedure must not only define the position of particles at initial time, but also associate to each particle a phase 
drawn from a uniform distribution. If the system is open, injected particles are given a phase when they enter 
the system. 

By definition of the shuffle updates, at each time step the particles are updated in order of increasing phases: 

Til < ^2 < •• • < Ti„, (3) 


where N is the number of particles present at a given time in the system. Note that with shuffle updates, if 
successive sites are occupied by particles with increasing phases, the set of particles can move as a whole one 
step forward within one single time step. 

Equivalently to the previous definition, in a continuous time picture, one can consider that each particle i is 
updated at time s + r^, where s is the discrete time. 

Two variants of shuffle updates have been considered in the literature. In the frozen shuffle update, phases 
are kept unchanged during the whole simulation. By contrast, in the random shuffle update, phases are drawn 
anew at the beginning of each time step. 

Shuffle updates have low fluctuations, in the sense that each particle is updated exactly once per time step. 
Still, for a given particle, the time between two updates can fluctuate between 0 and 2 for the random shuffle 
update, while it is always exactly one in the frozen shuffle update. Indeed, in the frozen shuffle update, a particle 
i is updated at times r^, + 1, + 2,.... 

When particles hop on a one-dimensional lattice with a constant hopping probability p, the model is equivalent 
to the so-called totally asymmetric simple exclusion process [1]. In this ID case, the fundamental diagram, i.e. 
the current J as a function of p, has already been calculated analytically for both updates in the deterministic 
case ip = 1). For the random shuffle update |241126j . the current with deterministic hopping reads 


J{P,P = 1) 


P 

p(i-p) 

2p-l 



1 


for p < 1/2, 
for p > 1/2. 


(4) 


In particular, this current has a maximum = 4 at p 
fundamental diagram reads 


JiP,P = 1 ) 


P 

2(1 - p) 


4. For the frozen shuffle update [231 HD) the 

for p < 2/3, , . 

for p > 2/3 '' 


with a maximum Jp^ = 2/3. This maximum corresponds to configurations where, after a short transient, 
particles form platoons. More precisely, an observer staying on a site and recording its occupation and the phases 
{xi} of the particles that occupy it at each time step would see a sequence of particles with increasing phases, 
then a hole, then another sequence of particles with increasing phases, and so on. A sequence of particles without 
a hole between them will be identified as a platoon. If we denote the average number of particles in a platoon by 
iz, then V particles will exit every v + 1 time steps, thus giving a maximal current For phases drawn from a 


^ See | 32| for a comparison of the queue shapes when Euclidian or Manhattan distances are used. 
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uniform distribution of values between 0 and f, v can be shown to be equal to 2 [23j . hence the maximal current 
value. 

We see that . Indeed, the repeated changes in the updating order of random shuffle update provides 

a supplementary source of blockings between particles, leading to less regular motion. 

For the frozen shuffle update, one must keep in mind that the flow may depend on the exact realization of the 
updating order, i.e. on the specific choice for the phases. This is the case for closed systems, or for the evacuation 
problem of this paper, or more generally for any system that involves a bounded number of particles whatever 
the simulation duration is. The phases are like a frozen disorder, and one needs to take ensemble averages to get 
the mean behavior [23]. 

By contrast, in open systems with incoming flows [28] . the time average would coincide with an ensemble 
average over the frozen disorder, i.e. over the phases, and a single simulation would be enough. 

2.4 The hybrid shuffle update 

Systems with random shuffle and frozen shuffle updates do not behave in the same way. For the first one, 
the renewing of the updating order introduces some randomness in the bulk particle dynamics while for the 
second, the disorder is injected at the boundaries and the bulk dynamics is deterministic (if the hopping rules 
are deterministic). 

It is also possible to combine both aspects by introducing a dynamics of the phases themselves. The choice 
of this dynamics is not unique. We will present one possible choice of hybrid shuffle updates in this sectioij^ 
This choice was guided by the request to have an update behaving similarly to the frozen shuffle update when 
the particle density is low, and closer to the random shuffle update inside congestions. We will discuss further 
our choice in section]^ but first we will define the hybrid shuffle update of this paper. 

We introduce the following rule. If a particle i attempts to hop from {x, y) to (x + a, y) (or (x, y + o)), and if 
the neighbours of the arrival cell (x + a, y ± 1) (or (x ± 1, y + o)) are both occupied by other particles, then the 
phase of the hopping particle Ti is redrawn from a uniform distribution on 0 < < 1 . 

The phase dynamics has no incidence on the motions of particles, i.e. the hops are always allowed when 
particles chose a neighbouring empty cell, according to the probability (i). Only the order of update will be 
modified, starting from the next time step. 

As particles on (0,0) are supposed here to quickly step out, and thus not to be in the congestion anymore, 
we do not redraw the phase when a particle hops to (1,0) even if (0,0) and (2,0) are occupied. 

Note that, in the one dimensional case, no renewal of phases occurs by definition. Therefore the hybrid shuffle 
update scheme is equivalent to the frozen one in one dimension, and we have 

Jl? = 4 ° = |. ( 6 ) 

As we mentioned, the shuffle update schemes presented in section |2.3| were already studied in detail in one 
dimensional systems (at least for the deterministic TASEP). Therefore one of our particular interests in this article 
is to investigate their behavior in a typical two dimensional system, i.e. an evacuation from a room though a 
small exit. In the next two sections, we shall explore the behavior of the evacuation model, and distinguish two 
dynamical regimes. 


3 Low-density regime 

At initial time t = 0, N particles are dropped randomly on the lattice of figure[^and their phases {xi} are drawn 
randomly from a uniform distribution. In this section and the next one, we show simulation results that we have 
performed with the room size = 51^ and various initial numbers N of particles, using either the random, 
frozen or hybrid shuffle update scheme. We will be interested in the total evacuation times T, i.e. the time when 
the last particle leaves the room, and the flux J of particles exiting the room, which is identical to the current 
between (0,1) and (0,0). 

Figure shows the total evacuation time T versus the initial number N of particles. When the number of 
particles increases, there is a smooth transition between two qualitatively different regimes, referred to as the 
low-density and high-density regimes. Let us first study the low-density regime. 

®In reference m we have introduced another hybrid shuffle update. The main difference is that here target sites are necessarily 
empty, hence the phase dynamics cannot be based on the occupancy of the target site as in reference m- The choice made in this 
article leads to higher flows in high density regime. 
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For low densities, the evacuation time is observed to be independent of the update scheme. Indeed, in this 
free flow regime, particles exit almost freely, and the velocity of non-interacting particles is the same for the three 
updates. 

The increase of the total evacuation time T with N is then simply due to the increase of the average of the 
distance between the farthest particle and the door. A particle starting on cell {x,y) is at a Manhattan distance 
d = |a:| -b y from the exit. Neglecting all interactions with other particles, and assuming k = oo, this particle will 
need a time t = d + 1 = |a;|-|-y-|-lto evacuate the room. Recalling that the shape of the room is square with 
cells {L = 21 + 1) and that the exit is located at the middle of one side, the number n{d) of cells located at a 
Manhattan distance d from the exit is found to be 


1 

f 2fi- 1 

for 

l<d<l, 

n{d) = < 

' 21 + 1 

for 

1 + 1 < d < 21 + 1, 

1 

61+ 4-2d 

for 

2l + 2<d<Sl + l 


Note that the n{d) takes a maximum value L when I + 1 < d < 21 + 1. Then we shall show that the average total 
evacuation time T is given by 


T = 1 + 



-1 3i-l-l 


d=l 






( 8 ) 


where 




is a binomial coefficient. In equation (8) the difference of binomial coefficients is the 


number of initial configurations where the farthest particle is at distance exactly d from the exit. One obtains 

the probability that the farthest distance is d, by dividing it by the total number of initial configurations f ^ 

and then the summation over d = 1,..., -|- 1 with multiplication by d gives the average farthest distance. The 

expression § of T versus N is plotted for L = 51 in figure]^ (a). We observe that the agreement is good when 
N is small, i.e. as long as particles (almost) do not interact. However, the true evacuation time becomes larger 
than predicted by (|^ as A^ increases, as a result of the collective effects occurring in the high-density regime. 

We now discuss in more detail the crossover towards this large N regime. Because of the geometry of the 
room, the number of particles arriving at the exit is not constant in time. A jam will form at the exit if the 
incoming flow becomes larger than the maximal current that the exit bottleneck can sustain. This limit 

value is called the capacity of the bottleneck, and is equ al to the outflow in the high-density regime, denoted by 
. The value of will be determined in subsection 4.1 For now it is enough to know that takes some 


numerical value that depends on the update scheme, but not on L or t. 

Let us first compute the incoming flow of particles. For the time being, we assume that particles move 
without interacting with each other. In the case fc = oo, the time t needed by a non-interacting particle to arrive 
at the exit site (0,0) is equal to the Manhattan distance d = \x\ + y between its initial position (x,y) and the 
exit. Then, for an initial uniform density N/V^ of particles, the inflow coming to the exit at time t is 


J“(t) = n{t) X N/L^, (9) 

where n{t) is given by Q. Between an initial increase and a final decrease, the function J*"(t) reaches a plateau 
value = N/L. 

When congestion occurs, i.e. one enters the high-density regime. The crossover thus occurs, 

when the number of particles is 

( 10 ) 

In the next section we shall focus on the high-density regime and determine the update-dependent value of 
This will allow us in particular to locate the crossover between the two regimes. 

For finite k, a similar transition is observed, see figure]^ (c,d). The form (|^ is no longer valid in the low- 
density regime for finite k. Instead of this we compared the genuine evacuation times of N particles (colored 
markers) with the single-particle problem {i.e. initially only one particle is put randomly in the room). The solid 
lines were obtained by averaging the maximum of evacuation times Ti of N independent simulation runs of the 
single-particle problem, i.e. formally 

E[max{Ti,.--,rAr}]. (11) 


4 High-density regime 

In this section we study the evacuation from the room in the high-density regime. 
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Figure 2: Total evacuation time T versus the number of particles N with random (red x), frozen (blue O) S'lid hybrid 
(green □) shuffle updates in linear (left) and logarithmic (right) scales. The room size is = 51^, and the parameter k is 
chosen as fc = oo (a,b) and fc = 3 (c,d). The plot markers were obtained by averaging over 10^ simulation runs. The black 
solid lines in (a,b) and (c,d) correspond to ^ and (111, respectively. For comparison, we drew dashed lines J = const, x t 
in (b) and (d). In (a) and (c), the arrows indicate the crossover points (101 between the low- and high-density phases, see 
the last paragraph of Subsection |4.1[ 


4.1 Outgoing currents {k = oo) 

From figure it can be seen that in the high-density regime the evacuation time grows linearly with the number of 
particles. Moreover, it can be checked in figure]^ that the particle current flowing out of the system (the outflow 
J]p) is constant in time except in the initial and final transients, so that the evacuation time is approximately 
given by T = N/J^, where J|P depends on the update schemes, i.e. X =R (random shuffle), F (random shuffle) 
and H (hybrid shuffle). It remains to compute the value of this current for these three updates, which is done 
for fc = oo in the following paragraphs. 

We start with the case of random shuffle update. The outflow will be determined by the dynamics around 
the exit. We therefore write an approximate master equation for the where the letters a and h stand for 

the occupations of cells (0,1) and (0,0), respectively, e.g. we write a = 1 or a = 0 depending on whether (0,1) 
is occupied or unoccupied. Since we expect that there is a queue around the exit, we suppose that (1,1), (0, 2) 
and (—1,1) are always occupied. Under this first order approximation and with the simple choice k = oo, cells 
(0,1) and (0, 0) are never simultaneously empty at integer times and we can write a closed equation for the three 
remaining probabilities as 

/ Plit) \ / 0 1 1/2 \ 

P{t +1) = MP{t), where P{t) = I P'i’(f) hM= \ 1/4 0 1/5 . (12) 

V Pllt) J \ 3/4 0 3/10 / 

Let us explain for instance how the first column of the matrix can be obtained. It corresponds to the transition 
probabilities from state J to states J, ° and J. Starting from state J, the particle in (0,1) will surely hop to (0, 0) 
during the time step. Consider the particles occupying sites (0,1), (1,1), (0,2) and (—1,1). If the particle on 
(0,1) draws the highest phase (probability 1/4), the three other particles will be blocked and the state at the 
next timestep will be ° whereas in the opposite case (probability 3/4) at least one particle from (1,1), (0,2) or 
(—1,1) will be able to hop towards (0,1) and the next state will be (. 

By calculating the stationary solution i.e. P = MP, we find a stationary current 

43 

J2D ^ po pi = _ = 0.60563 .... (13) 

This prediction is in good agreement with the simulations, see figure |^a). 
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Figure 3: Outflows as functions of time for random shuffle (a), frozen shuffle (b) and hybrid shuffle (c) updates 
for three values of A: = 1, 2 and 10 and for a room of size = 51^ with quarter-filling, N = — l)/4 = 650 

(the side L always has an odd number of sites). We are thus clearly in the high-density regime. The theoretical 
predictions for fc = oo are = 43/71 for random shuffle update and Jp^ = 1 for frozen shuffle. For comparison, 
we also drew the line = 2/3 for the hybrid case. 


Let us turn to the frozen shuffle update. For the 2D case, platoons are also formed similarly to the ID case: 
An observer staying on the exit site and recording its occupation and the phases {ri} of the particles that occupy 
it at each time step would see a sequence of particles with increasing phases, then a hole, then another sequence 
of particles with increasing phases, and so on. A sequence of particles without a hole between them will be 
identified as a platoon. Platoons are self-organized in a stronger sense than the ID case. The ordering of the 
particles is not conserved in 2D. And thus they become very large. Numerically we observe that the average size 
of these structures keeps growing with the number of particles in the room. This fact leads to the current 

= 1 (14) 


which is expected to become the exact value in the limit N ^ oo. 

In the high-density regime of the two-dimensional evacuation, the hybrid shuffle update is very similar to the 
random shuffle update because there is a queue near the exit. When a particle i hops, most of the neighbouring 
cells are occupied, so that its phase will often be redrawn. The outflow of the hybrid case is therefore similar to 
that of the random case, as shown in figure The remaining difference comes from the fact that, in the hybrid 
case, phases are often but not systematically redrawn. Still, as we shall see below, it is enough to strongly limit 


the 2D outflow as we wanted. While we were able to predict the flow precisely for the random case in (12), it 
is not so easy to do it for the hybrid shuffle update, because of the remaining ‘correlations’ of the phases. By 
simulation results, the two dimensional current for the hybrid case is measured as 




2D 


0.64. 


(15) 


The maximal one dimensional flow for the hybrid shuffle update scheme is equal to 2/3, as stated in (j^. 
The interactions between particles result from the dynamics of the phases {ri}. 

Now we go back to the problem of the crossover between low- and high-density regimes. By substituting 
the values of the outflows in the high-density regime (13), (14) and ( [T^ into (10), we estimate the transitions 
between the low- and high-density regimes to occur at N,, « 30.8,51 and 32.6 for the random, frozen and hybrid 
shuffle update schemes, respectively, for L = 51. These values agree with the numerical data, see figure(a). 


4.2 Dependency on k 

In the preceding subsection we studied the case k = oo. Here we consider the effect of varying k on the total 
evacuation time T. 

In figure]^ (a) we plotted the total evacuation time as a function of k for the three updates. The evacuation 
time is monotonous and decreases with increasing k. In this congested two-dimensional situation the hybrid 
shuffle case is closer to the random shuffle one, as expected. 

Now we apply an approximation, which is similar to what we did for k = oo: all the sites (x, y) with |a;|-|-|y| = 2 
are assumed to be occupied at any integer time. For the calculation of hopping probabilities, we also use an 
approximation that we take into account only the two choices, i.e. to stay at the same site and to move the 
neighbouring site closer to the exit. Then the probability of particles’ hops (0,1) —>■ (0,0), (—1,1) —(0,1), 
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Figure 4: (a) Total evacuation time per particle as a function of k and (b) its relative variation with k for the random (red 
x), frozen (blue O) s^nd hybrid (green □) shuffle update schemes with L = 51^ and N = (L^ — l)/4. For better visibility, 
we used logarithmic scales for k and T{k)/T(oo) — 1 in (a) and (b), respectively. In (a) we observe the agreement with 
the predicted evacuation time T/N = 1/J with equation (13l in the random shuffle case (dashed line). For comparison, 
we also drew the dotted line T/N = 1.5, which is obtained if we assume that the outflow is identical to the ID flow in the 
hybrid shuffle case. To obtain the data points we averaged over 10^ simulation runs. Due to the finiteness of simulation 
runs, fluctuations for large k are observed in the log scale plots. 


(0, 2) —>• (0,1) and (1,1) —>• (0,1) (denoted by po,Pi,P 2 and pa, respectively) are given as 


1 


^ — k 


r. — k 


PO = 


1 + e 


life, Pl=P3 = 


D— fc p—y/2k ’ 4“ 


(16) 


Now we extend the master equation (12) to general k. Since (0, 0) and (0,1) can be simultaneously empty, we 
add Pgit) to the probability vector: 


f nit) \ 
mt) 


nn 


P{t + 1) = MP{t), P{t) = 

V Plit) j 

For the random-shuffle case, the elements of the transition probability matrix are given as 

\ 


(17) 


M = 


S 3 

S 3 

0 

0 

0 

0 

3 +S 1 +S 2 + 3 S 3 

12 -^0 

i2+3si + 2^2 + 3 S 3 
60 

- S 3 

1 - S 3 

1 - To 

1 - 2 TO 

0 

0 

9-si-S2-3s,-i 

18—3si —2 s2—3s3 

12 Po 

60 


Po 


Po / 


(18) 


where Si = 4-92+ 93 , S 2 = 9 i 92 + 9 i 93 + 9293, S 3 = 919293 , qi = l—pi. Solving the balance equation P = MP, 

we finally find the outflow as 


Mk) = po + pi = 


_ 120po (1-S3)+Po(l-S3 )(9+si-52-9^3) _ 

120(1 —S3)+2po(42+3si+2s2—24 s3+2siS3+3s2S3 + 12s2)+Pq(9+si —S2— 9 S 3 ) 


(19) 


The total evacuation time is then estimated by 7 r(A:) = N/JYi{k). 

For the frozen and hybrid cases, the correlation among phases prevents us from obtaining M explicitly. In 
figure]^ however, we luckily observe that the plots Tii{k)/T}i{oo) — 1 and the plots T^{k)/T yi{oo) — 1 almost 
overlap. This means that the total evacuation time for the hybrid shuffle update with finite k is estimated as 


T}i{k) = rH(oo)J r(oo)/ Jpi(fc) by using (19). In the frozen case, when a particle of a platoon does not hop. 


another one may squeeze in and modify the structure of the platoons, hence a stronger effect of finite k values. 

Finally we discuss on the crossover points between the high- and low-density regimes. For finite fc, the inflow 
of the low-density regime at time t, J“(t), is not given by (|^. We rather determine it from single-particle 
simulations, using the approximation 

J“(t) « Af X P(Ti = t), (20) 

where P(Ti = t) is the probability of the evacuation time Ti in each simulation run of the single-particle problem. 
Then the crossover points are estimated from max J'”(t) = i.e. 


« J°'^7maxP(Ti = t) 


( 21 ) 


with the outflows measured in the three updates with finite k. For fc = 3, the crossover points are numerically 
estimated to be ~ 32.1,48.4 and 34.2, which provide good estimations, see figure [^(c). 
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5 Pedestrian evacuation 

5.1 Collective effects 

Some of the most impressive collective effects exhibited by pedestrian crowds are found in evacuations, see 
reference |34j for review. Countless human disasters have happened during mass events, where the motion of the 
crowd lead to stampedes [35113S] ■ 

Modeling is useful to predict the crowd behavior, but also as a tool to understand the link between individual 
behavior and collective effects. In particular, interactions between individuals play a key role in evacuation 
processes. 

Several types of interactions take place during evacuation processes. At very high density, the flow at a 
bottleneck can be completely stopped due to the formation of arches, in which the pedestrians are compacted 
and cannot move because of physical contact with their neighbours. 

At low densities, pedestrians interact without physical contact through so-called social interactions, in order 
to avoid collisions while minimizing the inconvenience of deviating from their optimal trajectory. A full charac¬ 
terization of these interactions is difficult, as they are non local, non isotropic, they include some anticipation, 
and they are non-additive when more than two pedestrians are involved. 

Another way to characterize the interactions is to observe their consequences on the collective behavior and 
to infer from their most relevant characteristics of interactions. 

5.1.1 Obstacle in front of the exit 

One interesting feature is the role of an obstacle suitably located near the exit of a room. Some experiments have 
been performed, which compare the outflows from the room, with and without obstacle. When pedestrians are in 
a panic-like condition, be. with strong physical contacts, the presence of a well located obstacle can surprisingly 
increase the outflow by about 30% [3S]. But for obvious security reasons, this effect has not been systematically 

confirmecQ 

In more “normal” conditions, the difference between outflows with and without obstacle are far less pro¬ 
nounced |35l [311 132] (the outflow increase is estimated around 4% in [ST], 6% in [32]) and it should still be 
confirmed whether this weak increase is also observed in real life conditions. However, even if the outflows are 
more or less equal, this is still surprising, as we could expect that an obstacle would decrease the outflow (and 
indeed this is the case for ill-located or ill-sized obstacles [32]) • 

5.1.2 Faster-is-slower effect 

Another effect that was conjectured from numerical observations is the so-called faster-is-slower effect, namely the 
fact that a stronger will of pedestrians to go out may lead to longer evacuation times. An experiment of aircraft 
evacuation showed a trend to have larger evacuation times in competitive trials compared to non-competitive 
ones, but the result was strongly geometry dependent and opposite effects could also be observed [33]. More 
recently, it was shown [inisi] that indeed a faster-is-slower effect could be obtained with humans, leading to an 
increase of evacuation times of typically 15 to 20% under competitive egress, and related to physical contacts and 
the formation of clogs. A faster-is-slower effect can also be observed in some other systems that exhibit clogging, 
namely granular matter or sheep |41] . By contrast, as ants do not clog but rather keep their density always more 
or less constant even under stress [32], no significant faster-is-slower effect was observed for ants [33], except in 
the trivial limit where the stress source has an impact on the physiology of the ants, and thus on their stepping 
behavior |42j . In this paper, we rather focus on normal evacuation conditions of pedestrians, be. with limited 
physical contacts, for which no faster-is-slower effect was ever evidenced. It is an open question to characterize 
the interactions between pedestrians in these more ordinary situations. 

5.1.3 ID and 2D outflows 

An experimental observation that may seem less spectacular but has strong implications for modeling, is that, 
if one compares the outflow of a single line and of two-dimensional flows, one finds that 2D outflows are similar 
or even sometimes smaller than ID outflow^ Indeed, in m, two experiments are performed and the reported 

^ A similar effect was observed for ants m and for granular matter m, but the systems are too different to extrapolate to 
pedestrians. 

^This remark applies for disordered random 2D flows. If pedestrians are pre-ordered into several incoming lines, higher flows can 
be achieved EH ED. but this situation is quite special and we shall not consider it here. Still, it is an open question to know how 
to account for this phenomenon in simulations. Due to the discrete lattice, straight motions in directions not parallel to the lattice 
have to be decomposed into hops of different orientation. Besides, the exclusion rule is too crude to account for the subtle zipper 
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decrease of the 2D outflow compared to the ID outflow is respectively 4% and 14%. 

The fact that 2D outflows are not larger than ID outflows is quite surprising. In fact, in two-dimensional 
flows, more pedestrians arrive at the exit. This overfeeding effect should lead to higher flow^ Actually this is 
what we observe when we compare the ID and 2D outflows obtained in sections |2.3| and : Both for the frozen 
and random shuffle updates, the 2D outflow is larger than the one-dimensional one. This makes it more clear 
now why we have introduced the hybrid shuffle update defined in section 2.4 One sees from figure]^ (c) that 
for this update, the 2D and ID outflows (both measured for the deterministic case fc = oo to compare equivalent 
models) are similar. 

The fact that 2D outflows are not larger than ID outflows indicates that there are some interactions between 
pedestrians beyond the hard-core ones, that compensate the overfeeding. As these interactions decrease the flow, 
they are often assimilated with friction. One must be careful that this notion of friction can refer to different 
types of interactions (with physical contacts or without in ‘normal’ flows), with different consequences on the 
collective behavior as we explained above. The choice of the friction implemented in the simulations implicitly 
assumes whether physical contacts are considered or not and it is not clear yet whether a single model can account 
for all these different regimes. 


5.2 Updates 

5.2.1 Parallel update 

Most cellular automata models for pedestrians use parallel update. This update has the particularity that 
conflicts occur: when two pedestrians have the same target site, the exclusion rule requires to select which one 
will actually hop. The need to include a conflict resolution procedure is in general considered as a drawback in 
simulations. However, in the case of pedestrian flows, it was suggested that these conflicts could have a physical 
relevance [HI [301 mi- Indeed, they occur preferentially in these converging areas where friction should occur. 
An easy way to tune the friction level is to allow with a certain probability none of the pedestrians involved to 
move. This clearly reduces the overall flow when the number of conflicts increases. 

5.2.2 Shuffle updates and stepping 

One interesting feature of the frozen shuffle update is the possibility to interpret the phases as the phases in the 
stepping cycle of pedestrians. This allows in particular, for ID deterministic flows, to map the cellular automaton 
dynamics in the free flow phase onto a continuous space and time dynamics [28] , reproducing the regular motion 
of free pedestrians. This also gives a physical interpretation to the priority given to one pedestrian when two 
pedestrians meet. The time step of the update, which plays the role of a reaction time, is then directly connected 
to the stepping period - actually it is known that pedestrians cannot turn/change velocity at any points of the 
stepping cycle but only at specific moments [45) . In the frame of this interpretation of the phases, it becomes 
natural to modify the phases in the jammed phase, as done by hybrid shuffle updates. Indeed there is no reason 
why a pedestrian should keep its stepping phase while being blocked by its predecessor. 

It should be noted that cellular automata are rather mesoscopic models that do not aim to reproduce the 
precise microscopic dynamics but simply to give the correct average flows and densities at specific key locations. 
Still, it would be interesting to see if features coming from the stepping behavior of pedestrians could be retained 
and give interesting properties to the models, closer to what happens at the scale of individuals. 

The hybrid shuffle update presented in this paper illustrates that indeed some interactions leading to flow 
reduction in the congested phase can be implemented through a partial renewal of the phases. The choice we 
have made is not unique (see |33j for another one). One difficulty is to trigger friction preferentially in convergent 
flows - and not in all high density situation^ Still, we have not explored all the possibilities opened by hybrid 
updates. While in this paper the phase dynamics consists simply in a partial redrawing of the phases from a 
uniform distribution, one could also consider some evolution rules for the phases where the new phase value 
would depend for example on its previous value, on the phases of its neighbours, or possibly on the phase of the 
last visitor of the target cell. 

5.3 Fundamental diagrams 

Flow-density relations have been measured in several experiments. However a consensus is difficult to reach, 
first, because the fundamental diagram depends on the geometry of the system, and second, because the precise 

effects that are favored with pre-ordered lines at the bottleneck m- 

^ Indeed, as we said in the previous footnote, this increase is what happens with pre-ordered lines |31lI32| . 

^ For the parallel update, conflict occurrence automatically coincides with convergence of the flow. 
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Figure 5: Typical queue shapes for random shuffle (a), frozen shuffle (b) and hybrid shuffle (c) updates with 
parameters L = 51^, N = (L^ — l)/4 and k = 5. 


method used to measure velocity and density may have a strong influence on the results naiis]. However, now 
that many data have been collected, one can nevertheless make a few observations. 

A first set of experiments deals with one-dimensional flows, where pedestrians walk in a line without passing 
each other sziiiiiin]. If one would crudely approximate the fundamental diagram in this one-dimensional case 
by a triangular shape, the density for which the maximal flow should be clearly smaller than half the maximal 
density m- This is not the case for the one-dimensional TASEP, for which the density of maximum flow is 
p = 1/2 for random shuffle and parallel updates, and p = 2/3 for frozen shuffle update. 

It is known that one way to shift the maximum of the phase diagram towards smaller densities is to allow 
particles to jump more than one cell at each time step. Indeed this was shown for the Nagel-Schreckenberg model 
for car traffic 0131, but also for the floor field model with parallel update inilll]. 

In two-dimensional flows, the variety of possible geometries (corridors, bottlenecks, intersections, etc) makes 
it more difficult to conclude. In some measurements the density of maximal flow was found to be shifter 

towards higher densities than in ID, while in some others not. Besides, pedestrians can have different walking 
directions, a feature which of course impacts the fundamental diagram. 

In the simulations, the level of discretization m, or the details of the hopping rules [S3], can also modify 
the shape of the fundamental diagram. Thus it is not possible to validate a model only on the basis of its 
update scheme, though the update scheme has an influence on the properties of the model. The whole set of 
characteristics of the model should be taken into account. 

5.4 Queue shapes 

We end this discussion by a few observations on the queue shapes. Though having a realistic queue shape in 
front of the bottleneck is not necessarily required as long as outflows are correct, it would be more satisfactory 
to approach the shape of real queues. In figure 0 we show the positions of particles at different moments of 
time, as obtained with the model of this paper for A: = 5. The queues look like half-circles for the three updates 
considered. This is what is usually observed in simulations. In experiments, shapes like droplets have rather 
been observed (see figure 4 in [54] 1. This difference of queue shape between simulations and experiments was 
already mentioned in [34] . We have shown in reference [33] a simulation snapshot exhibiting a queue with droplet 
shape, as a result of some modifications of the rules, especially adding diagonal hopping. Indeed, it was shown 
in [55] how the queue shape can vary strongly (including very unrealistic shapes) depending on the allowed hops 
and their associated rates. Another way to decrease the density on each side of the exit is by penalizing sharp 
turns [31] . 
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Random 

Frozen 

Hybrid 


1/2 

2/3 

2/3 

j2D 

43/71 

1 

0.64 


Table 1: Values of the maximal currents for the random shuffle, frozen shuffle and hybrid shuffle updates in one 
and two dimensions (and for k = oo). All the values of are exact in the limit of an infinite system. The 
value of for frozen shuffle update is also expected to be exact in the limit of an infinite system. For random 
shuffle update, = |y was led by an approximation, and the 2D value for hybrid shuffle update was measured 
numerically. 

6 Conclusion 

In this article, we have investigated the properties of shuffle updates in two-dimensions through the study of a 
simple evacuation problem. A particular attention has been given to evacuation times and outflows, which we 
have studied both numerically and analytically. 

In this study, we have considered a floor field model with only a static field. In general floor field models, a 
dynamic floor field is also considered, which allows to transform some long range interactions between pedestrians 
into purely local ones [16]. Indeed our aim was not to design a full pedestrian model but to explore how 
interactions between pedestrians could be modified through the update scheme itself. More generally, it is 
an important issue to understand how to model properly inter-pedestrian interactions, as they determine in 
particular the outflow of evacuation problems. 

We have used a unified picture to define the various shuffle updates, by associating a phase Ti to each particle. 
Then updates differ only by the dynamics of these phases. For the random shuffle update, all the phases are 
redrawn in every time step. For the frozen one, phases are initially (t = 0) given to particles, and we keep them 
till the end of evacuation. For the hybrid one, we partially redraw some of the phases. The last one interpolates 
between the frozen at low density and the random in congested configurations. Using the hybrid shuffle update 
the phases {ri} keep their interpretation as trackers of the walking cycle, and the dynamics of the phases in the 
high density phase simply corresponds to a loss of memory of the stepping cycle phase when pedestrians are 
forced to slow down or stop. 

It is clear that many updating procedures that interpolate between frozen shuffle and random shuffle exist, 
and that the choice we have presented in this article is only a particular instance of them, that we have chosen to 
illustrate how model properties can be modified by including some phase dynamics. We found that this test case 
has similar ID and 2D outflows. Table shows the flows for the shuffle updates, summarizing subsections 
12.41 and l4.1l 

Other features could be included, in particular via the phase dynamics. For example, in reference [55] 
measurements of the time gap between two consecutive pedestrians are presented. In our model, phases allow to 
have an interpretation of the dynamics in continuous time, and hence it is possible to define time gaps. Attempts 
could be made to monitor the distribution of this observable, e.g. by drawing the phase of a blocked particle 
in correlation with the phase of the blocking particle. 

The idea of deriving an updating order from the stepping dynamics was even pushed further in the so-called 
optimal steps model [531551155], for which not only the phase but also the duration of the stepping cycle and the 
length of the steps vary from one individual to another one. Though in the case of a constant step length the 
optimal steps model can also be applied to cellular automata models, it is more generally defined in a continuous 
space. It would be interesting to study the influence of this update on outflows too. 

Improvement of pedestrian models would require a better understanding of flow-density relations. Beyond 
the fact that a better consensus should be reached on which fundamental diagram the models should verify, one 
must be aware that fundamental diagrams are supposed to be obtained in stationary regimes. In practice, in 
experiments, one realizes short-term and localized averages. In the extreme case where instantaneous individuals 
measurements are realized, one observes quite different flow-velocity relations. Indeed, the fluctuations around 
the mean behavior can be quite strong |49j . It is an open question to determine whether real pedestrian flows - 
including the important case of evacuations - can be fully described by implementing stationary laws, or whether 
they could be dominated by transients, which can give quite different flow values, not necessarily well reproduced 
by models calibrated in the stationary regime. Another issue to improve cellular automata models would be to 
study more systematically the isotropy of their properties m- 


2.3 
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